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which is justified in a large Nf approximation. The back reaction on 
the Bose fields due to fermion field fluctuations is calculated using a 
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for the response of fermions to sphaleron transitions, the renormal- 
izability of the effective equations of motion and non-perturbative 
dynamics in the framework of nonequilibrium quantum field theory. 
The long-time behaviour of the system is discussed and we speculate 
about applications to finite density calculations. 



*email: aarts@phys.uu.nl 
■f email: jsmit@phys.uva.nl 



1 



1 Introduction 



The time evolution of quantum fields in and out of equilibrium plays an impor- 
tant role in the early universe and in heavy ion collisions. Detailed questions 
concerning baryogenesis, inflation, the chiral phase transition and time evolution 
during phase transitions in general, require a non-perturbative analysis, prefer- 
ably directly in real (as opposed to imaginary) time. Since exact answers are 
usually not available, a considerable amount of work has been devoted to finding 
suitable approximation methods. 

One approach which has received a lot of attention recently is based on Hartree 
or large N approximations |2|, [3]], and has been applied to pair production in 
strong electric fields [|J, (inflationary) cosmology || || 0]j the formation of chiral 
condensates in heavy ion collisions || and nonequilibrium symmetry breaking 
0. In this approximation method, the fields are written as the sum of 'classical' 
or 'background' fields (which are then called order parameter, or condensate^) 
and quantum fluctuations around them. A set of coupled equations for the back- 
ground fields and the fluctuations can be derived from first principles in several 
ways from the Schwinger-Keldysh formalism for real-time (non-) equilibrium field 
theory, and they have the form of semiclassical equations of motion. A nice prop- 
erty of the equations is that they contain only the usual quantum divergences 
which can be renormalized with the usual counterterms. 

However, in numerical work the emphasis in the literature has been mainly 
on situations where the classical field is homogeneous. Although this is sufficient 
for some situations (inflationary scenarios, pair creation in strong homogeneous 
electric fields), this is certainly not the case for sphaleron transitions or realistic 
treatments of heavy ion collisions. 

A useful toy model for obtaining a better understanding of the dynamics that 
plays a role in various scenarios for baryogenesis (for recent reviews, see, e.g., 
[|IIJ) has been the abelian Higgs model in 1 + 1 dimensions [Fl, The model 



has many properties in common with the electroweak sector of the Standard 
Model, in particular it has a non-trivial vacuum structure which is labeled with 
the Chern-Simons number, and sphaleron barriers separating the vacua. Fur- 
thermore, extending it with axially coupled fermions results into a model with 
anomalous fermion number non-conservation, similar to baryon number violation 
in the Standard Model. 

In this paper we study this abelian Higgs model with axial fermions, using 
the nonequilibrium method mentioned above. We treat the bosonic fields as 
classical background fields, and the fermions as fluctuations, which is justified 
for Nf fermions in the limit of large Nf. The bosonic fields are not restricted 
to be homogeneous, allowing for genuine sphaleron transitions. The equations of 
motion for the bosonic fields have the form of classical field equations coupled 



1 We shall use the terms classical field, background field and order parameter interchangeably. 
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to the fermion current or force. The fermions are treated with a mode function 
expansion, in which the mode functions obey a Dirac-like equation in the presence 
of inhomogeneous background fields. 

We formulate the model on a lattice in space and time. This allows for a 
simple and stable time integration of the coupled partial differential equations 
for finite lattice spacing. Furthermore, gauge invariance can be implemented and 
exact symmetries of the lattice action give rise to exact conservation laws at finite 
lattice spacing. 

In the literature, to our knowledge fermion fields have been treated numer- 
ically only in the case of homogeneous mean fields {|, |7[] (pair production, pre- 
heating in inflation) , while fermions in the presence of inhomogeneous mean fields 
have sofar been investigated analytically only (see e.g. [13]). 



In the remainder of this Introduction, we give the outline of the paper. In 
section 0, we introduce and summarize some basic properties of the abelian Higgs 
model with axial fermions in the continuum. In particular, we recall that the 
axial gauge symmetry can be transformed into a vector one, which is convenient 
because it allows us to extend known methods for fermions on the lattice in 
imaginary-time to real-time. This transformation leads to an unusual Majorana- 
type Yukawa coupling, hence, we go from a formulation with Dirac spinors to one 
with Majorana spinors. The next two sections are rather technical. In section 
|3j we present the model on the lattice. Fermion fields on a lattice lead to the 



fermion doubling phenomenon [14 , 15] , 16], and we describe in this section how 



we deal with this. Symmetries, currents and useful observables are discussed in 
section |j. The large Nf effective action is briefly derived in section |5| and the 
resulting equations of motion are given in section [| The effective equations of 
motion are solved numerically and we present results in sections |7|-|T0|. In section 
[7] we use handmade time dependent Bose field configurations that go from one 
vacuum to another via a sphaleron transition. These external fields are used to 
test the properties of the fermion system. Issues related to renormalization are 
discussed in section l| In section P, we present numerical results for the full non- 
perturbative dynamics, including the fermion back reaction on the Bose fields. 
We discuss the long-time behaviour of the system in section |1(] and also mention 
a possible application to finite density. We end with a summary and outlook. 

Conventions and technicalities are summarized in several appendices. Ap- 
pendix |A] contains our conventions and in appendix [B| we give the initial mode 
functions used for the fermions. Finally, in appendix [C] we show how the rela- 
tion between spectral flow and the anomalous conservation law works out in our 
lattice model. 
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2 Continuum model 

We consider the abelian Higgs model in 1 + 1 dimensions with axially coupled 
fermions. The action consists of a gauge, Higgs and fermion part and is given, in 
the continuum, by 

S = S G + S H + S F , (2.1) 

with 

S G = - J d 2 x^F^F^, (2.2) 

Sh = - j ' d 2 x[(D,4>yD^ + V(4>)}, V^) = X(4>*4>-v 2 /2) 2 , (2.3) 

S F = - J d 2 x + iqfasW + G$(<f>*P L + <f)P R )^} . (2.4) 

Here F^ u = d^A u - d v A^ = <9 M - %A^ and P R>L = (1 ± 7s)/2. The model 
is defined on a circle with circumference L in the spatial direction. The bosonic 
(fermionic) fields obey (anti)periodic boundary conditions. Our conventions are 
summarized in appendix |X[ 

The model is invariant under local axial gauge transformations 

V>(x) -> e^ hs ^(x), (j){x) -> e-*<*ty(a:), -> - d M f (z). 

Gauge invariance implies that the fermions have one half of the charge of the 
scalar field, q — |. The conserved (axial) current, = 0, associated to the 
global gauge symmetry, is 

f = -J 5 M + ft, 3s = ni>-fl^ f h = i(D»<py<t> - icj>*D^cj>. 
Furthermore there is a global vector symmetry with a classically conserved current 

corresponding to fermion number. 

The global vector symmetry is anomalous, and the anomaly equation reads 

9fj,jf = —qto P - 

The topological charge density can be written as the divergence of the Chern- 
Simons current 

gtop = 4^F Mi/ = d,C\ C" = ^A v . 
Fermion number violation is expressed as 

Q f (Jt)-Q f {0) = -(C(t)-C{0)), (2.5) 
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relating the total fermion number and the Chern-Simons number 

Q f (t) = J dx x f f {x), C{t) = J dx l C°{x) = --^ J dx 1 A 1 {x). (2.6) 

It is possible to transform the model to one with a vector gauge symmetry, 
by performing charge conjugation only on the right-handed fields This is 

convenient, because we want to treat the model on a lattice.^ The transformation 
is explicitly 

?P R =(^ R C) T , 4> R = -(C^' R ) T , ^ = V4, i>L = $ L - (2-7) 

C denotes the charge conjugation matrix. In these primed variables, the fermion 
part of the action is 



Sp = — / d x 



{21 



In the Yukawa mass term only the antisymmetric part survives. The gauge- 
fermion part of the action, with G = 0, has now the form of the Schwinger model, 
or QED in 1 + 1 dimensions, with q = \ charged fermions. The Yukawa term, 
coupling the fermion field to the scalar field, has turned into a Majorana-like 
mass term. 

In terms of the primed fermion field, the gauge symmetry is vector-like 

- e-^y(x), 4>{x) - e-^V(z), A,(x) - A^x) - d^(x), 
with the conserved current 

f = jf + f h , ff = *#V¥ = -it 

The global symmetry that involves only the fermion field is now axial, with a 
classically conserved current 

f ^ e -^^', ft = $VW = ~ff (2-9) 
The anomaly equation reads 

Q' 5 (t) - Q' 5 (0) = C(t) - C(0), Q' 5 (t) = [ dx l j' 5 °(x). (2.10) 



In the rest of this paper we will only work with the transformed action. Therefore 
we drop the primes from the fermion fields and charges in the following sections. 

In the remainder of this section, we summarize some properties of the model. 
We work in the temporal gauge A = throughout the paper. The bosonic 



2 A similar transformation can also be performed in the electroweak SU(2) Higgs model in 
3 + 1 dimensions. 
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vacuum can be labeled with the Chern-Simons number C. The local energy 
minima are characterized by C being integer valued and \<f>\ = v/y2. The vacua 
are separated by finite energy sphaleron barriers, where C is half-integer. The 
height of the energy barrier is given by the sphaleron energy E sp ^ = ^\/2Xv 3 
We will show the vacuum structure in more detail, also in the presence of 
fermions, when we use handmade sphaleron transitions in section [7|. 

The dimensions of the fields and coupling constants is as follows: <p an d v are 
dimensionless, ip has a mass dimension of 1/2, A M ,e and G have dimension 1, 
and A has dimension 2. 

Elementary powercounting shows that besides the partition function, only the 
one loop self energy of the scalar field is superficially divergent. This divergence 
can be canceled with the appropiate choice for the bare scalar mass (or v 2 ). 

Finally, the tree level masses, due to classical Higgs expectation value, of 
respectively the scalar, gauge and fermion fields, are 

G 2 v 2 

m 2 H = 2\v 2 , rn 2 Ah = e 2 v 2 , m 2 F = — - — , (tree level). 

At the one loop level, the fermion field not only affects the scalar self energy, it 
contributes also to the mass of the gauge boson 

Q 2 C 2 

m 2 A f = , (one loop fermion). 

7T 

This concludes the introduction of the model. In the next section, we dis- 
cretize it on a lattice in space and time. 



3 Lattice formulation 

In section |] and beyond, we derive and study a set of effective equations of 
motion, that involve the Bose fields and the back reaction of the fermion field. 
This study is done numerically, and involves a discretization. A convenient way 
to discretize the equations is by use of a lattice formulation in space and real 
time. There are good reasons to use a lattice: exact symmetries of the lattice 
action give rise to exactly conserved quantities (i.e. for finite lattice spacing), 
the lattice distance is a gauge invariant ultraviolet cutoff, the resulting time 
integration algorithm is simple and stable. For a general introduction to lattice 



gauge theories in euclidean spacetime, we refer to ||16|| . Before we write down the 
effective equations, we first discuss in this section the lattice formulation and in 
the next section exact symmetries and conserved currents. 

The model is defined on a lattice in space and (real) time, x = (x°, x 1 ). 
When there is no confusion we also use a; = t, x 1 = x. The lattice spacing in 
direction \i is denoted by a M = (a , a\ = a), and there are iV spatial lattice points 
(L = aN). As indicated in the previous section, from now on we only work with 
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the transformed action with a vector-like gauge symmetry (i.e. ( [2.1| - [2.3| , |2.8| ) in 
the continuum), hence we drop the primes from the fermion fields and currents. 

To put the Bose fields on the lattice is straightforward. We replace the (co- 
variant) derivative by the lattice (covariant) derivative as 



[U^x)f{x + a^)-f(x)} 



dj(x) = - [f(x + a,) - f(x)} , Dpfix) = 

d 'J( x ) = — lf( x ) ~ f( x - ^)] > D 'nf( x ) = — [f( x ) - u »( x - ^)f( x - ^)] 



a, 



a, 



denotes the unit vector in the \x direction, and we use the notation 

The gauge and scalar part of the action keep the same form as in the continuum 

Sg = -VAV. (3-1) 



(3.2) 



X,fJ, 



We use the non-compact formulation for the gauge field. 

To put fermions on a lattice is less straightforward. The naive way leads 
to so-called fermion doublers in space and time 0, [15| . Instead of one physical 
fermion particle, the lattice action describes, close enough to the continuum limit, 
2 2 physical fermions. One effect of the doublers is to cancel the anomaly, when 
the naive continuum definition of Q§ is used In the euclidean formalism, 

the doublers can be removed by adding a higher derivative term to the action: 
the Wilson term, proportional to a parameter r [TJJ]. Another possibility is to 
interpret the doublers as physical particles, i.e. the lattice action describes more 
than one flavour. 



19] 



This is the idea behind the staggered fermion method 
In this paper, we add a Wilson term in space, proportional to n, to suppress 
the space doublers. We interpret the doublers in time as physical particles. Hence 
the lattice action describes two flavours. In |2(J we discuss in detail fermions on 
a lattice in real-time, and we construct the Hilbert space and the transfer matrix, 
associated to the lattice fermion path integral, which clarifies further the physical 
interpretation of the lattice system. 

The fermion part of the lattice action is given by 



X 

The Dirac operator is 



i/>(p + W)i/) + -Gty VCty - -G#Cf 
2 2 



£7*A«> A. = 5 to* + ^i) 



(3.3) 
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which is explicitly (we recall that the fermions have charge q = |) 

D, xy = gi- [6 x+ ^, y U^x) - 8 x , y+ ^{y)} , Ufa) = e~^ A ^\ 

Note that = — Z? M , and the (anti) symmetric parts are given by = ImD^, D A 
ReD^. The Wilson term (in space and time, with parameters = (ro, is 

= [~25 x , y + 6 x+iXa ^Ufa) + 6 x , y+ ^U$(y)] , 

with W /t = W, and the (anti)symmetric parts W s = ReW, W A = ImW. As 
indicated above, we use r = throughout this paper. The numerical results 
presented further on, are obtained using r\ = 1.Q 

Because of the Majorana coupling between the scalar and the fermion field, 
it is convenient to use a Majorana formulation with real gamma matrices and 
a four-component real Majorana field It is constructed from the real and 
imaginary parts of i/j, 

^ = ^(* 1 -i* 2 ), ^ = -^=(^ + z^), * = ( ^ ) ■ 

We also introduce Pauli matrices, written as pi^s, that act on the 1,2 components 
of*. 

In terms of * and the pi matrices, the fermion part of the action can be 
written as 



X 



A + P2 p S + W S + P2 W A + -^=(0lP3 + faPl) 



Here we defined the covariant derivative and the Wilson term on a Majorana field 
as 

i.e. 



In the original model of the previous section, the Wilson term corresponds to a Majorana- 
Wilson term 

Sf,w = -J2 [V^PlW^ - i>P R CW$ T ] . 
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qp2 is the Majorana charge matrix. For the Yukawa part we used the explicit 
representation C = ft, 

and we introduced a matrix field notation 

$ = ^|(0lP3 + 02Pl)- 

The lattice field equation (in the temporal gauge) reads 

^i{do + & )^ = n D ^, (3.4) 
where the Majorana Dirac hamiltonian is given by (a 1 = — 7°7 1 ) 

H D = -ia 1 V 1 +f3(W 1 + G$). (3.5) 

It is hermitian and antisymmetric (Hd^ = T~Ld = —7~Ld T )- 

The lattice fermion doublers in space are removed by the Wilson term in 
space. We now describe how we deal with fermion doubling due to the 'naive' 
time discretization (ro = 0). Reverting to ordinary Dirac notation, the free 
fermion lattice propagator is given by 

r, t \ m F P - Hid' 1 sin ap + i^aT, 1 sin a uj 

S{P,U) = 2,-2-2 =T—2 • ( 3 - 6 ) 

m F + a z sin ap — a sin clquj 

The fermion mass term mp p = mp + m p is the sum of the physical continuum 
mass term mp and a term due to the higher derivative in space, the Wilson 
'mass term', m p = a -1 Ti(l — cos ap). This propagator takes the usual continuum 
form for a, a — > 0. However, the propagator has also a pole near a^uj = ±7r, 
corresponding to the rapid time variation (— 1)*/°°. As shown in ||15|| , such poles 
have to be interpreted as particles, the doublers. Their physical frequency p° is 
related to u by a u = n + a Q p°, such that we recover the continuum propagator 
in the limit a, a — > 0, up to the unitary transformation ift-f 5 = —ji, 



7r 0n I m F - + ry u p L 



S(p, — + P u )^7i Y-— 2 — 7i- (3-7) 

a V mp+p 2 -p z J 

Hence, the lattice fermion field ip describes two physical fields in the continuum 
limit, which we label with flavour indices u, d, and denote with Xu and Xd, 

^(x,t) - X u(x,t) - (-lf ao llXd (x,t). (3.8) 

In Majorana notation, x — {Xi — iX%) / v2, this takes the form 

V(x, t) -> X u {x, t) - (-lf ao llP2 X d (x, t). (3.9) 
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Substitution in the lattice action displays its physical content in the continuum 
limit 

S F -> - J d 2 x X -X T i5 [fDp + G$r 3 ] X. (3.10) 

Here we combined X u ^ into a flavour doublet X on which Pauli matrices 71,2,3 
act. Cross terms oc (— l)*/ a ° do not contribute in the continuum limit. A useful 
relation is 

7iP2^d(+G) 7 iP2 = -n D (-G). (3.11) 

The Majorana Yukawa coupling happens to break the SU(2) flavour symmetry 
of the continuum limit of the G = theory. Experience with euclidean lattice 
fermions suggests that it could be avoided, if so desired, and that we can obtain a 
wide class of target continuum theories in a more sophisticated staggered fermion 
formulation. 

It is convenient to formulate the back reaction of the fermions on the Bose 
fields (in sections ^] and beyond) in terms of initial expectation values of fermion 
operators. This calls for an extension of ( |3.9|) to an operator equation. Since 
the continuum fields X u ,Xd are supposed to be slowly varying in lattice units, 
they are somewhat spread out over the lattice. A minimal spreading in time is 
over two time slices, which is also used in derivations of the quantum mechanical 
Hilbert space from the lattice path integral pi]] . Labeling pairs of time slices by 
an integer k, we may write without 'overcounting', for even t/a = 2k, 

= X u (x,t) -jip 2 X d (x,t), 
V(x,t + a ) = X u (x,t)+^p 2 X d (x,t), (3.12) 

This relation suggests that it is also useful to combine ^(x, t) and ^/(x,t + a ) in 
a doublet, 

B<^) = 4 ( •&tV° ) Y (313) 



V2\ *fat) 
so that ( p,12| ) can be written concisely as 

5 = RX, (3.14) 

with 

R=^=(] 7lP2 ), Rrf = h RR T = n. (3.15) 

Eqs. ( |3.12j ) (or equivalently ( (3.14|) ) can now be viewed as relations between 
operators. In |2^| we will give more details on the Hilbert space aspects of fermion 
doubling in real time. For now we assume that the physical Majorana fields X 
correspond to hermitian operators satisfying the anticommutation relations 

{X(x,t),X T (y,t)} = 6 x , y . (3.16) 
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Note that as operators, the \I/'s are not hermitian because of the imaginary p^. 
Typically, it is the Grassmann variables ip{x, t+ao) and ip + (x, t) which correspond 
to hermitian conjugate operators in Hilbert space ]3T|, |20| . 

We can expand the fields X in eigenspinors of the Dirac hamiltonian in flavour 
space, 

H m = -ia 1 V 1 +P(W 1 + T 3 G$), (3.17) 

and identify the coefficients in the expansion with creation and annihilation oper- 
ators. For example, let U au (x), U a d(x) be a complete orthonormal set of positive 
energy^ eigenspinors of Hdr at t = 0. Here a labels the set (which typically 
contains momentum), and u, d indicates the flavour. We then expand 



X U (X, 0) = J2 KuUau(x) + blJjl 



IX 



(3.18) 



and similar for d. The creation and annihilation operators obey the usual anti- 
commutation relations 

{baf,b a 'f'} = S aa 'Sff, f = U,d. (3.19) 

We can now express the lattice field \1/ at time slices t/ao = 0, 1 in terms of the 
eigenspinors and the operators, if we combine (|3.18| ) with Q3.12| ) at k — 0. This 
gives 

9(X,0) = Y^[ b -uUau(x)+bl a U* au (x) 

a 

-badl 1 p2Uad{x) ~ b^ 1 p 2 U* ad {x 

*(x,Oo) = [ b <*uUau(x) + b]JUx) 

a 

+b ad l l p2Uad{x) + b ] a( fi l p 2 U* d (x] 

This relation can be extended to arbitrary times by writing 
tffot) = KuU au (x,t) + bl u U* au (x,t) - b ad U ad {x,t) + bi d U* d (x,t)] , (3.20) 



where U au (x,t),U ad (x,t) are solutions of the lattice Dirac equation ( |3.4j ), with 
the initial conditions 



U au (x,0) 
U ad (x,0) 



U au (x,a ) = U au (x), 
-U ad (x,a ) = "i l p2U a d{x). 



(3.21) 



^The negative energy solutions are given by U* 
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We shall refer to the U(x, i)'s as the mode functions. We show in |2IJ that solving 
the Dirac equation for the mode functions U is equivalent to time evolution in 
Hilbert space by means of the transfer operator. Note that (|3.21|) shows that 
the ordinary modes start smoothly in time, while the doubler modes start like 

The mode functions remain orthonormal under the discrete time evolution in 
the appropriate inner product. In terms of 

V af (k) = ( Ua u af ~[ t ) 0) ) > t = 2 ^o, f = u,d, (3.22) 

where we suppress the space dependence, the time evolution can be written as 
V af (k + 1) = M(k)V af (k), 

' ' x " ~ \ -i2a n D2k +i 1 

where the subscript k indicates the time dependence of the Bose fields in the 
Majorana Dirac hamiltonian. The matrix M is unitary in 'the r± inner prod- 
uct', MWiM = T\ and it follows that the norm of the V's is conserved in this 
inner product. The orthonormality relations follow from the initial conditions: 

(r 3 )///. Similarly the anticommutation relations of the 
lattice Majorana field are preserved: {E(x,k),E T (y,k)} = Ti5 Xjy , which follows 
directly from (|3.14|-|3TTED. These conservation relations have the same meaning 



as the Wronskian condition ||-[[|. They imply that the basic anticommutation 
relations, ( |3.16 ) and Q3.1M) , remain valid also for times larger than zero. 



To summarize this section, we have formulated the model of section g on a 
spacetime lattice. A large part has been devoted to the physical interpretation 
of the lattice fermion field. It represents two physical fermion fields in the con- 
tinuum. In relation with the dynamics to be treated in sections [| and beyond, 
the important result of this section is the mode function expansion of the lattice 
field (|3.20|) in terms of time independent operators (|3.19|) and spinor mode func- 
tions. These mode functions are solutions of the Dirac equation Q3.4D with initial 
conditions (|3.21| ). 



4 Symmetries, currents and observables 

We discuss several exact symmetries of the lattice action, since these lead to 
exactly conserved currents for finite lattice spacing. Furthermore, we construct 
other observables that are related to symmetries of the continuum theory. 

The local gauge symmetry under which the action is invariant, acts on the 
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fields as 

A^x) -> A^x)-d^{x), U^x) -> e-^f/^e^+H 

We indicated how the fields transform for the different ways in which they are 
written. 

The global gauge symmetry leads to the conserved current j^(x) = j^(x) + 
jf{x), with the scalar part 

j£(x) = i{D»4>{x)y<j>{x) - t(f)*(x)D^(f)(x) } (4.1) 

and the fermion part (in Dirac and in Majorana notation) 

f f {x) = tq[^(x + a^P^U q J(x)^(x)-^(x)P^(x)^(x + a^)] (4.2) 
= % - [V T (x + ap)pP$q P &${x)V{x) - y T (x){3P^qp 2 U«(x)y(x + o M )] . 

We use the notation, familiar from euclidean lattice fermions, P± = (r M ± 7 M )/2 
(recall that r = 0). It is easy to check what this current represents in terms 
of the flavour fields X U ,X^ using the relations ( |3.12| ). For example, the charge 
density at even t/a is 

jf(x,t) = ^[$ T {x,t + a )qp 2 y{x,t) + y T {x,t)qp 2 y{x,t + a )} 

= \ [ X u( x , t)qp 2 X u (x, t) + Xj(x, t)qp 2 X d (x, t)] , 

which is precisely the charge density due to two Majorana fermions with charge 
q, as expected. 

There is another global symmetry, due to the fact that the model is defined 
on a lattice. It is inspired by the staggered fermion formalism. It is 

V{x,t) ^exp [iu^ 1 {-l) t/a °] #(z,t). 

The associated charge density is 

( — l)t/ao 

j°{x, t) = ^ [^ T (x, t)^{x, t + ao)- ® T {x, t + a )^(x, t)] . 

The physical interpretation becomes clear in terms of the flavour fields. At even 
t/ao, 

3a(x, t) = \ [X%(x, t)p 2 X d (x, t) + Xj(x, t)p 2 X u (x, t)] 
= -X T (x,t)p 2 r 1 X(x,t), 

13 



where the last line is written compactly in flavour space. It is the first component 
of the flavour triplet current. 

Other observables, which are not directly related to a symmetry of the lattice 
action, are constructed with the continuum limit as a guidance. The energy 
density of the fermions is found to be 

H f (x, t) = \ [® T (x, t)H D (x, t)V(x, t) + (t^t + oq)] . (4.3) 

In terms of the flavour fields at even t/ao, using again ( |3.12| ), this is the expected 
energy density 

H f {x,t) = ^X T {x,t)H Dfl {x,t)X{x,t), 

where Hdr was defined in ( |3.17|) , and we used (|3.11|) . The anomalous axial 
vector current which satisfies the anomaly equation is the flavour singlet current, 
i.e. where the contributions from both flavours add, 

j° 5 (x,t) = ^X T (x,t)^P2X(x,t) 

= ~ [X*(x, t) l5 p 2 X u {x, t) + Xj(x, t) lb p 2 X d {x, t)] . 

It is straightforward to verify for even t/ao, that in terms of the lattice field the 
axial charge density is obtained by taking 

f 5 (x,t) = J [^ T (x,t)j 5 p 2 ^(x,t) + (t^t + a )] . 

We will divide, in the rest of the paper, the axial charge (density) by a factor of 
two to keep the anomaly equation, given in section ||: Q$(t)— Qs(0) = C(t)—C(0). 

The current associated with the global symmetry ip —>■ exp (— iu^) ip of the 
continuum action, is given by 

jS&W = \ [4>(x)^ l5 U«(x)iP(x + a,) + 4,(x + a,)^ l5 Uf(x)ij(x)} , 

= % - [^ T (x)(3Yl5P2^(x)^(x + a M ) + V T (x + a^lsU^(x)^(x)] . 

The continuum interpretation of this current is not that it represents the anoma- 
lous axial vector current, but rather the third component of the flavour triplet 
axial vector current. For example, the charge density at even t/ao is 

j° 5 (x, t) = \ [^ T (^ t + a )j 5 p 2 y{x, t) + V T {x, t)j 5 p 2 V{x, t + a 
I 

= 2 [ X « ( x ' t h5p2X u (x, t) - Xj(x, t)j 5 p 2 X d (x, t)] 
= ^X T (x,t) l5 p 2 r 3 X(x,t). 
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Indeed, this density represent the difference between the two axial charge densities 
instead of the sum. This is known in the naive fermion formalism as well . The 



current is not exactly conserved because the global symmetry is broken explicitly 
on the lattice by the Wilson term in space. 

After this lengthy discussion of the model and its lattice version, we come to 
the effective action and equations of motion in the following section. 



5 Large Nf effective action 

In this section we derive the effective action by solving the path integral to leading 
order in the large Nf approximation. Since this can already be found in great 
detail in the literature |2|, [|, we only show the essential steps. The effective 
equations of motion follow directly from the effective action. 

To avoid a notational jam, it is convenient to denote the Bose fields collectively 
with ip = ((p, <f>*, Afj). We also denote the corresponding sources as J = (J*, J, J M ). 
We use a compact notation, such that e.g. J ■ ip = ^ J(x)<p(x), where the sum 
is over all indices, including space-time. 

The partition function we are interested in is 

Z [J] = Tr pTc exp iJ ■ ip. (5.1) 

Tq indicates that the exponent is time ordered along the Keldysh contour in the 
complex time plane, which is shown in figure [I]. We take the initial density matrix 
p such that the fermion field has zero mean value. Matrix elements of the initial 
density matrix are written as 

((fa, ^ a \p\Pb, *&) = p(<Pa, * a, <Pb, 

The corresponding path integral expression reads 

Z[J]= J Vip a Vip b V^ a V^ b p{ip a ,^ a ;ip b ,^ b ) J V<pV*exp(iS + U-(p), (5.2) 

with S = Sb + Sp, Sb = Sq + Sh- The c again denotes that the whole expres- 
sion is defined along the Keldysh contour. The initial values of the fields and 
correlations are determined by the density matrix. 

Let's for a moment consider the part of the partition function that involves 
the fermion field. Formally, we can perform the integration over the fermion field, 
and we denote the result as 

^log/%] =log J V* a V%p{(p a ,* a ;<p b ,* b ) Jv* expiS F [(p]. (5.3) 

D[ip\ depends on the Bose fields ip along the contour, including the initial Bose 
fields ifaifb- Note that if the density matrix is quadratic in D represents the 
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Figure 1: Keldysh contour in the complex time plane. 



resulting determinant. We would like to remark that it is not necessary to have 
a gaussian density matrix for \l/ to perform the integration. Furthermore, we do 
not need to specify any property of the density matrix regarding its Bose field 
dependence. 

To arrive at a large Nf expansion, we duplicate the field ^> Nf times, and 
demand that all fields ty l are initialized with the same density matrix, i.e. 

p{ip a ^ a -ip b ^ b ) -+ ]Jp(ip a ,¥ a ;ip b ,¥ b ), i = l,...,N f . 

i 

If we then integrate out the fermion fields, we find 

Z[J) = J Vcp a V(p b j Vcp exp % {^S B - i^Nf \ogD + J ■ cp^j . 

To proceed, it is common use to rescale the parameters, and the Bose fields and 
sources such that the exponent in the path integral becomes proportional to Nf. 
In our model, this rescaling is explicitly 

= (j>'y/N~f, V = V'y/N}, J = J'y/Wf, 

e = e'/y/Nj, G = G'/y/Nf, A = \'/N f . 

Note that this rescaling implies that the original, unsealed theory is weakly cou- 
pled when Nf is large, and that the initial conditions for the unsealed <fr and A^/e 
have to be such that they are of order y/Nf. 

In terms of the primed variables, the partition function is (up to normaliza- 
tion) 



Z[J'] = J WaW b J W exp iNf (s' B - i\ log D' + f- cp'^j . 



(5.4) 



which suggests a saddle point expansion for large Nf. The field configurations 
that determine the stationary phase, (p s , are found by extremizing the exponent 
of 0. 

A Legendre transformation from the generating functional log Z depending on 
the sources, to the effective action depending on the so-called classical or mean 
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fields, 



Z[J'\ = exp iN f W[f], <p = r[fl = W[J'] — J'-ip, 

yields the effective action, which is to leading order in large Nf, 

Y[ip\ = S B [ip\- % -\ogD[(p\. 
The effective equations of motion follow from 

d£M = _ f 

dip 

The mean fields <p are defined on the complete Keldysh contour. In general, 
the sources J are different on the upper and lower part of the contour (which 
is needed to be able to construct correlation functions of arbitrary observables) . 
This implies that also the mean fields are different. However, having derived the 
equations of motion and expressions for the observables, we can set the upper and 
lower sources equal, which implies that the upper and lower mean fields are also 
the same. The initial conditions are in principle to be derived from the density 
matrix. This can be analyzed conveniently if it has a gaussian form. To leading 
order p s = (p [|] . 

Differentiating log.D at t > t{ gives the back reaction of the fermion field in 
the background of the Bose fields, at time t. We use a bracket notation to denote 
the result 

1 d _ dS F 

2d<p(x,t) og ~ l{ d<p(x,ty' 

The initial conditions for the fermion field are also determined by ( |5.3|) . However, 
we find that is more transparent to use the operator formalism. We discuss the 
fermion initial conditions in the following section. 

To conclude this section, we would like to remark that in this approximation 
scheme the fermion fields are integrated out exactly. The large Nf approximation 
is only used to justify a saddle point expansion of the resulting effective bosonic 
action. The leading order result in this expansion is that the Bose fields are 
treated classically, and the fermions show up in a one loop contribution, but with 
the exact propagator. In the following, we drop the bars from the mean fields. 

6 Effective equations of motion 

We are now in a position to give the effective equations of motion. Using the 
results from the previous section, we find the bosonic field equations (in the 
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temporal gauge) 

d' d A 1 {x,t) = e 2 [jl(x,t) + (jj(x,t)}] , (6.1) 

d , 1 d A 1 {x,t) = -e 2 [j° h {x,t) + (j%x,t))} , (6.2) 

d' Q d (f){x,t) = D[D l( f>(x,t) - 2\[\0(x,t)\ 2 - 4/2] <j){x,t) + G(F(x,t))(6.3) 

The scalar and fermion contribution to the current, j£ and j*i respectively, are 
given by (|4.1| ) and ( [4.2|) . The force on the Higgs field due to the fermions is 
defined by 

t] = ha^j) = 4* T( * ! mpi + lp ^^ t] - (6 - 4) 

Eq. ( |6.2| ) is Gauss' law. It implies (on a periodic lattice) that the total charge is 



zero, 



Q(t) = Q h (t) + Q f (t) = OhM + (j%x,t))) = 0. (6.5) 



In ( |6.3| ), we used a subscript b to indicate that v 2 B is a bare parameter. We'll 
discuss renormalization in section || 

Following for the case of homogeneous classical fields, we calculate the 

back reaction of the fermions on the Bose fields using a mode function expansion, 
introduced in section ||. We expand the field \1> in a complete set of eigenspinors 
of the initial Dirac hamiltonian, with time-independent creation and annihilation 
operators as coefficients (i.e. ( |3.20| )) 

*(x,t) = [bauU au {x,t) + bl M U* u {x,t)-b ad UU^t)+bi d U* d {x,t)\ . (6.6) 

a 

Expectation values of the operators determine the initial quantum state of the 
fermion field. Two typical choices are a vacuum and a thermal state. These are 
given by the following expectation values 

(b a fb a >f) = n af 5 aa >5f f > {b a fb\, f ,) = (1 - n af )5 aa >5f f >, f = u,d. 

Here n a f is the Fermi-Dirac distribution 

n af = [exp(E af /T f - m ) + l]" 1 , (6.7) 

with E a f the eigenvalue of the Dirac hamiltonian for eigenspinor U a f, and Tf^ n 
the initial temperature of the fermions. We will often encounter 

(Kf,bl f }) = 1 - 2n af = a af . (6.8) 

All fermion observables can now be expressed in terms of cr a f and the mode 
functions. Regarding the numerical results that are presented in this paper, we 
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restrict ourselves to initial fermions in vacuum, i.e. 7/i n = 0, for which n a f 
vanishes of course. 

The system is closed by giving the equation of motion for these mode func- 
tions. As explained in detail in section || both U au and U a d are solutions of the 
lattice Dirac equation ( |3.4j ), 

^i(do + d' )U af (x,t) =n D [A 1 (x,t),4>(x,t)}U a f(x,t), f = u,d, (6.9) 

with initial conditions ( |3.21|) . The number of partial differential equations that 
has to be solved is quite large. Besides the Bose fields, there are AN four- 
component complex mode functions, which make the numerical calculation scale 
as iV 3 (for fixed Qq/o). 

To be explicit, we list the fermion field expectation values in terms of the 
mode functions. To avoid too lengthy formulas, we combine the two flavour 
mode functions in one, 

TT I U aU 

a ~{ u ad 

and use for the expectation values ( |6.8| ) the matrix 

a m , 



U -cr ad 

We stress however that the mode functions are still to be solved with the original 
Dirac hamiltonian (i.e. not the one in flavour space). 

In the equations of motion we need the charge and current density, and the 
force the fermions exert on the scalar field. These are given by 



(f f (x, t)) = -j^Y] Kfo t)qp 2 a a U a (x, t + a ) + U^(x, t + a )qp 2 cr a U a (x, t 

a 

(jj(x,t)) = i -J2[Ui{x,t)(3P 1 _qp 2 cr a U q 1 {x,t)U a {x + a,t) 

a 

-Ul(x + a,t)/3Plqp 2 a a Uf\x,t)U a (x,t)}, 
(F{x,t)) = - l -J2[Ui{x,t)P{p 1 + tp 3 )a a U a {x,t)}. 

a 

Two observables we are interested in, are the anomalous axial charge and the 
conserved flavour charge. Their densities are 



(j° 5 (x, £)) = *)75P 2 a Q £/ Q (x, t) + (t -> t + oq)] , 

a 

(ja(x,t)) = ^[Ul(x,t + a )'j 1 a a U a (x,t)-Ul(x,t)'j 1 a a U a (x,t + a )]. 
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Finally, the energy density of the fermions is 

E f (x, t) = -- [Ut(x, t)H D (x, t)a a U a (x, t) + {t^t + a )] . 

a 

All these formulas are valid for a complete set of mode functions that are 
initially eigenspinors of an arbitrary Dirac hamiltonian, i.e. a Dirac hamiltonian 
with arbitrary Bose fields as background. In the numerical treatment of the 
equations we present in the following sections, we use as initial background fields 
for the Dirac hamiltonian a vacuum configuration of Bose fields, i.e. A% — 0, (j) = 
where vr is the renormalized vacuum expectation value. In terms of the 
density matrix, this means that we take no initial correlations between the Bose 
and the fermion fields, p = ps <8> Pf, and that pf is quadratic in In Appendix 
|B| we calculate the eigenvalues and -spinors in such a configuration, with the 
following results. 

The energy eigenvalues can be summarized as 

Epn = \J s l + m OT> 

with 

s p = a -1 sin pa, m pv = m p + rjmp, r\ = ±. 

Here p is discrete and takes N values, given by (|B.2| ). The total fermion mass is 
the sum or difference between the mass due to the Wilson term, m p = a~ Y r i(l — 
cos pa), and the mass due to the Higgs vacuum expectation value, mp = Gvr/ a/2. 
The collective label is now explicit, a = (p, rj), and has 2N values. For these initial 
conditions, we find that the charge and current densities vanish initially, 

(f f (x,0)) = (j}(x,0)) = (f B (x,0)) = (j°M)> = 0. 

The initial fermion force on the scalar field is 

(F(x, 0)) = (F(x, oo)> = ± vvn^r 1 - (6-10) 
The initial energy density 

E f (x, 0) = -- Y a Pv E pv 

PV 

contains the zero temperature divergent part (a pv — > 1). We renormalize this by 
subtracting the bare fermion energy density, 

Ef(x,t) = E f (x,t) - Ef, Ef = ~J2 E w ( 611 ) 

PV 
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The total energy of the system E tot is the sum of the energy in the fermions and 
in the classical Bose fields, 



E tot {t) 



E b (t) + E f (t), E b (t) = J2 le 2 E 2 (x,t) + \n(x,t)\ 2 + V pot (t), (6.12) 



X 



with 



e 2 E(x,t) = doAi(x,t), ir(x,t) = do4>(x,t), 



V P ot(t) = [\Di4>(x,t)\ 2 + X (|0(x,t)| 2 - 




X 



Let us conclude this section with the following remarks. When the initial state 
of the system is homogeneous, the system will remain homogeneous during the 
time evolution. The space dependence of the mode functions is given by plane 
waves for all times, and there is no mode mixing. The equations of motion reduce 
to ordinary differential equations. In the case of inhomogeneous initial conditions, 
the subsequent space dependence is determined by the equations. There is mode 
mixing, and the label p looses its interpretation for t > 0. Since the initial state of 
the fermions we use here, is choosen to be homogeneous, the inhomogeneity in the 
initial conditions has to be put into the Bose fields. This should be done in such 
a way that the total Higgs charge is zero, Qh(0) = (since Q/(0) = 0) and that it 
does not violate the local Gauss' law (|6.2j ). In this one space dimensional model, 
this can be done in a straightforward manner. Another possibility to introduce 
inhomogeneity would be to start not with a plane wave basis for the fermions, 
but with wave packets. This could be used to study scattering between fermions. 
This would correspond to an initial density matrix that is not quadratic in the 
fermion field. 

Before diving into the full dynamics, we first analyse the response of the 
fermions to certain external Bose field configurations. 

7 Handmade sphaleron transitions 

In the previous section, we listed the effective equations of motion and the ini- 
tial conditions for the fermion mode functions. Once the initial Bose fields are 
specified, the equations can be solved (numerically). In order to gain experience 
with the equations for the mode functions, we treat in this section the Bose fields 
as external fields. We focus on sphaleron transitions: the Bose fields start in a 
vacuum configuration, characterized by a certain integer Chern-Simons number 
C and scalar field winding number n, go through a sphaleron, and end up in a 
vacuum where C and n have changed by one unit. This analysis is closely related 
to the spectral flow analysis, which we summarize in Appendix 0. However, when 
considering spectral flow, usually the Dirac hamiltonian is diagonalized for every 
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external Bose field configuration p3, Oj. Here we use the real-time approach 



and solve the equations of motion for the mode functions in the presence of the 
Bose fields. Such a calculation can be done analytically in the case of QED with 



a homogeneous gauge field p3| . We do not have to discuss renormalization of 
the effective equations of motion at this stage, since the Bose fields are external, 
and the divergence is only present in the fermion back reaction. The divergence 



in the fermion energy is renormalized as in fl6.11|) . We stress again that the main 



purpose of this section is to test the equations for the mode functions. 

The Bose configurations we use are the following. A vacuum configuration 
with Higgs winding number n is given by 

<p n = -j= exp (2nin^j = A=V n . 

Note that on the lattice n is only defined modulo N, because Nx/L = x/a is 
integer. We use the gauge invariant definition of the winding number |24) on a 
periodic lattice, 

<f>{x,t) = pix^y 9 ^, 
n & = ^^([^ + M)-#(M)-^i(M)L + ^i(M)), 

X 

[fl = /(mod27r)e(-7r,7r]. 
The corresponding gauge field configuration is 

Al = f, C=-n. 

The sphaleron configuration is given by 

<p s = --^ tanh i^/\j2v{x - L/2)J exp (in- J , A x = --, C= -. 

We go from one vacuum to another through the sphaleron in a time to, with 
the sphaleron being halfway, to sets the time scale. Since we want to go through 
several vacua after each other, the configurations are taken to be (t' = t/to) 

4>{t') = /(t')0n„ + (1 - f(t'))<P s V nsJ A^t') = ^g(t'), 

where the winding number n v changes from n v to n v + 1 at t' being half-integer, 
and n s from n s to n s + 1 at t' being integer. At t' = 0, n v = n s = 0. The profile 
functions /, g that determine the time-dependence obey the following conditions: 
/ equals 1 at integer t', and at half-integer t', and g equals i! at (half-) integer 
t' . Furthermore, we require the time derivatives of /, g to vanish at t' being 
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(half-) integer This means that the configuration is 'stationary' both on top 
of the sphaleron and in the groundstate. Note that these configurations are 
inhomogeneous in general. 

The parameters we use in this section are 

N = 32, a /a = 0.05, eL = 3.2, et = 2, , . 

w 2 =4, A/e 2 = 0.25, G/e = 0.1, [ ' 

except when indicated otherwise. Let us shortly discuss the lattice spacing de- 
pendence. All the results we present are for a fixed physical size eL. Hence, 
increasing the number of lattice points N reduces in fact the lattice spacing a. 
Furthermore, the time scale eto = 2 we take here is sufficiently large such that 
all physical time scales are large compared to lattice scales, but also not so large 
that the adiabatic approximation (eto — > °°) tells all. 

In figure ||, we show a characterization of the Bose fields configuration, i.e. 
the Chern-Simons number, the winding number, and the potential energy in the 
Bose fields, normalized with the sphaleron energy. The potential energy in the 
Bose fields is the well-known periodic potential. 

In figure |3| we show the response of the fermions, i.e. the change in the axial 
charge, when going through the series of configurations of figure |[ The axial 
charge follows the Chern-Simons number, which changes from to 4. This is the 
anomaly equation. There are high frequencies oscillations present in Q$, which 
are more clear in figure |], where we show a blow up of figure [3] for two values 
of the lattice spacing in time, a$/a = 0.1, 0.05. These high frequency oscillations 
reflect that the lattice Dirac equation is a three term difference equation in time. 
Indeed, in the case of the free Dirac equation we find the solutions 

t) = e ipx - l ^ puJ , uj = ±E p , u = -t E p . 

The first solution represents the normal physical particles, the second solution is 
the sum of a high frequency part (oc l/a Q ) and a low energy part representing the 
doubler. The high frequency part gives a (— l)*/ a ° behaviour with an amplitude 
that vanishes as clq . 

Like the axial charge, the fermion energy has rapid oscillations, with an ampli- 
tude that decreases for decreasing ao/a. The energy transferred to the fermions 
and the sum of the fermion and the potential energy in the Bose fields is pre- 
sented in figure || We use here a very small ao/a = 0.005 to suppress the rapid 
oscillations. The fermions lift the degeneracy of the bosonic groundstate, and 
introduce local minima. The picture is symmetric for C — ► —C (i! — ► — t'), we 
have no CP violation in our model. 

We would like to stress that taking the doublers into account leads to a con- 
sistent time evolution. One can imagine to start with only the mode functions 

5 To be explicit, we use f(t') = cos 2 nt',g(t') = t' - {^y 1 sin47rf'. 
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for the normal particles. Let us recall that these are initialized as smoothly as 
possible, U au (x, 0) = U au (x, a ). Right after t = 0, only the normal particles are 
present. However, due to the non-linearity of the equations of motion, high fre- 
quencies will be generated, and in order to make sense out of these, the doubler 
interpretation has to be introduced after all. Hence, it is better to include them 
from the start. In terms of evolution in Hilbert space, including the doublers 
explicitly leads to unitary time evolution in the Hilbert space that is associated 
with the lattice path integral. 

Another test is the periodicity of Q5 on the lattice. As explained in Appendix 
for vanishing Yukawa coupling, if C — N(—N) the axial charge takes its 
maximal (minimal) value, and if C = ±2N, the total axial charge is zero again. 
Although this phenomenom is a lattice artefact, which is due to the finite number 
of states, it still serves as a check on the fermion system. We show the outcome 
of this test in figure ^| for iV = 16,32 and G/e = 0. Plotted are Q5/N versus 
t'/N, because the curves should then fall on top of each other for large N. We 
see this to be the case already for these small N values. The maximum Q 5 also 
agrees well with the value 2/tt ~ 0.637 for N — ► 00 (cf. (|C.4| )). Of course, this 
maximum is way out of the region where there is continuum behaviour, i.e. where 
the anomaly equation is satisfied (Q5 ~ C ~ t'), which appears to be \Q$\ smaller 
than N/4. 

We proceed with an investigation of the effect of increasing the Yukawa cou- 
pling. In figure |7|, we show the behaviour of the axial charge for three values 
of the Yukawa coupling, G/e = 0,0.1,0.5, up to a large Chern-Simons number 
C = N. In the continuum regime \Q$\ < N/4 = 8, the anomaly equation is 
(approximately) satisfied for finite G as well, but the deviation of Q5 from C 
becomes more and more important for larger G. This can be further investigated 
by taking a fixed finite coupling G/e = 0.1, and decreasing the lattice spacing 
(increasing the number of lattice points). The result is shown in figure |8|, for 
N = 16, 32, 64. We conclude that for a finite Yukawa coupling lattice artefacts 
are more important, and that for fixed Yukawa coupling, an increase in the num- 
ber of lattice points shows convergence to the continuum limit. It is also known 
from previous (lattice) studies (see e.g. p2[ ) that a finite Yukawa coupling has a 
relatively large influence on the relation between the anomalous charge and the 
Chern-Simons number. 

Summarizing the results obtained in this section, we found, using a controlled 
series of external Bose fields configurations, that the system behaves for zero and 
small Yukawa couplings as expected. For bigger Yukawa couplings, the lattice 
artefacts are substantially larger. With this in mind, we can continue and study 
the complete dynamics. 
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8 Renormalization 



We go on with the full dynamics, and include the back reaction. The full quantum 
theory contains ultraviolet divergences, which are regulated by the lattice cutoff. 
In this section, we discuss the renormalization and show numerical results to 
demonstrate its applicability. 

Let us first make some remarks about the role of the parameter v 2 . By 
rescaling according to 

<j> = v<t>', e = e'/v, G = G'/v A = X'/v 2 , (8.1) 

v 2 appears in front of the bosonic part Sr of the effective action, with the 
fermionic part — ^ilogD unchanged. Hence, if we also scale the bosonic initial 
conditions accordingly, the influence of the fermions is scaled by a factor 1/v 2 , 
and the fermion back reaction gets reduced at larger v 2 . This can also be seen 
directly from the effective equations of motion.^ This is of course not surprising 
since 1/v 2 is indeed the dimensionless semiclassical expansion parameter. 

As indicated in section |2], only the scalar self energy is divergent in the full 
quantum theory. The divergence is renormalized by the appropriate choice of v 2 B . 
Therefore, let's take a closer look at the equation for the scalar field ( |6.3| ). Because 
we treat the Bose fields as classical and only the fermions are integrated out, we 
expect to find only the fermion one-loop divergence of the Higgs self energy. To 
see this explicitly, we consider for simplicity a ground state configuration for the 
Bose fields, A\ — 0, = vr/\/2, where Vr is the renormalized, effective Higgs 
expectation value. Substituting this in ( |6.3Q , gives the gap equation 



Kv 2 R -v 2 B )^ = G(F)(v R ), (8.2) 
where we indicated that the mode functions are initialized with v — Vr. Using 



( 16.101) for (F)(vr), we see that it depends on Gvr through m m and E m . This 
makes ( |8.2j ) a highly non-linear equation in G and Vr . We can, however, perform 
a lowest order approximation, and expand the right hand side in G, to find (at 
zero initial temperature a pr) f = 1) 



G»jil r 1 / m 2 

,f|H = Tfi?^ 1 "^ 

where E p q is independent of Gvr. The sum contains the expected logarithmic 
divergence, which we indicated as logiV. Using this approximation, the gap 



6 To be precise, the rescaling is with a finite v = Vr, this introduces a rescaled bare v' B 
vb/vr. 
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equation flS.2| ) has two solutions 



, G 2 1 ^ 1 / ml 



The log divergence in the second solution is canceled with v\. 

In practice, we fix v\ and then find (for certain N) the corresponding bare 
parameter v B , simply from (|8.2|) , 

= — (F)(v R ), (8.3) 

A V R 

with (F) given by ( |6.10| ). We demonstrate this with numerical results. We use 
the following parameters 

iV = 16,32,48, oo/a = 0.1, eL = 3.2, . , 

w| = 11.15 A/e 2 = 0.25, G/e = 0.5, [ } 

and as initial Bose fields configuration, we take 

<p(x, 0) = — ^L, do(fi(x, 0) = 3e cos — — , 
v 2 £ 

0) = 0, d Ai(x, 0) = Gauss' law + e 2 . (8.5) 

We start with a slowly varying wave for do4>- The space dependence of OqAi is de- 
termined by Gauss' law ( |6.2|) (it is zero in this simple case), only the homogeneous 
part can be freely chosen (it is the 'conjugate variable' for the Chern-Simons num- 
ber degree of freedom). Note that for this type of initial conditions, the initial 
energy of the whole system is entirely in Bose fields with long wave lengths. 

In figure |9|, we show the Chern-Simons number as a function of time for three 
values of the lattice spacing (N = 16, 32, 48). C undergoes a damped oscillation. 
The bare expectation value is kept fixed, v 2 B = 10, and not changed according to 
( |8.3|) . The three curves do not agree, and the period is lattice spacing dependent. 



In figure ITU], we perform the same calculations, but now with a fixed renormalized 
expectation value, v R = 11.15. The bare expectation value that enters in the 
dynamics is determined according to (|8.3j ), and is N dependent. The value of 
vr is chosen in such a way that for N = 48, v 2 B = 10, such that the curves can 
be compared with the N = 48 curve in figure |^. In this case, we see that the 
three lines fall on top of each other, and there is convergence towards the physical 
continuum limit. The line that shows the largest discrepancy is for largest lattice 
spacing (N — 16). 

To show that the observed convergence is also true for other observables, we 
repeat the analysis for \(fi\ 2 = L~ x J2 X l < ^( x ; ^)| 2 - I* 1 figure PI, v 2 B is kept fixed, and 
the curves show again lattice spacing dependence. In figure O, we keep again 
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v\ fixed, and determine v\ from (|8.3|). Again this leads to converging physical 
results, with the largest discrepancy for smallest N, N = 16. 

We conclude that there is no convergence when the bare parameters are not 
adapted in the proper way, i.e. according to ( |3.3|) . A correct renormalization leads 
to convergence towards the physical continuum limit, already for fairly small N. 



9 Non-perturbative dynamics 

In the previous section, the total energy of the system was, in some sense, not 
very large. The Chern-Simons number simply oscillates around zero, with an 
amplitude which is not larger than 0.2. There are no sphaleron transitions, C 
does not get larger than 0.5. In this section we consider systems with more energy. 

To achieve this, we use the same type of initial Bose fields configurations as 
in the previous section, but with two wave lengths and complex amplitudes 

Vfi 27TX AlTTX 

(p(x, 0) = —=, d (p(x, 0) = 6e(l + i) cos — — h 2e(l + 2i) cos — — , 
V 2 L L 

Ai(x,0) = 0.1/L, a Ai(x,0) = Gauss' law + e 2 . (9.1) 

The parameters are 

N = 32, ao/a = 0.05, eL = 3.2, , , 

4 = 8(4), A/e 2 = 0.25, G/e = 0.0(0.1). { > 

We begin with the case of zero Yukawa coupling. The fermions start in a 
vacuum state, with zero energy. All the initial energy is in the Bose fields. We 
show the energy in the Bose fields Ef,, in the fermion field Ef and the total 
energy E tot in figure [13|, for fairly larger times than in the previous section, 
< et < 3000. The total energy is conserved up to small oscillations, very 
similar to the result from a leapfrog algorithm. The initial state is clearly a 
nonequilibrium situation, and there is energy transfer from the Bose fields to the 
fermion degrees of freedom. Presumably there is fermion particle production, as 
in the nonequilibrium homogeneous models studied before [f§]-|§. 

In figure [14] we show the Chern-Simons number and the axial charge. Initially, 
there is a lot of energy present in the Bose fields. The Chern-Simons number is 
wildly moving. The Bose system does not feel the sphaleron energy barriers. 
However, there is continuous energy transfer towards the fermions, and plateaus 
become visible. The energy transfer continues and the time the system spends 
in one 'vacuum' becomes longer. follows C, and it is not always possible to 
distinguish Q§ from C. This is actually a nice property, because it shows that the 
anomaly equation is correctly reproduced. Small differences are due to the finite 
lattice spacing. In figure |T5] we show the Chern-Simons number and the Higgs 
winding number. The winding number is anticorrelated with C when sphaleron 
transitions dominate, as expected from section |7J. 
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Let us this stage discuss the stability of the numerical algorithm. The total 
charge Q = Qf + Qh is conserved better than O(10 -10 ) (actually, for zero Yukawa 
coupling, both Qf and Qh are separately conserved). The first component of the 
vector triplet charge is conserved better than C(10~ 16 ). Locally, Gauss' law 
is respected better than 0(1CT 10 ). 

For finite Yukawa coupling, we found in section [7| with the help of handmade 
sphaleron transitions, that the axial charge is more sensitive to the lattice dis- 
cretization. In particular, for the typical parameters used in this section and with 
N = 32, the axial charge suffers quite a lot from discretization effects. This indeed 
shows up in the dynamics. In figure [H] we show C and Q$ for v 2 R = 4, G/e = 0.1 
and et between and 500. The anomaly equation is reasonably satisfied only for 
a finite time, in contrast to the G/e = example. 

10 Finite temperature and density 

Let us summarize what we have done so far. In all the results we have shown, 
the initial Bose fields are long wave length fields. The fermion contribution is 
calculated starting from a vacuum state, i.e. (|6.7|) with T/ iin = 0. It is straight- 
forward to start with a non-zero fermion temperature, and study e.g. the effect 
of (Landau) damping on the Bose fields due to a thermal bath of fermions. 

However, in this section we would like to comment briefly on the possibility 
of applying the methods used in this paper to a situation where also the Bose 
fields are thermal. Such a system would describe a thermal plasma consisting of 
classical Bose fields and quantum fermion fields. This would be interesting e.g. 
from the perspective of baryogenesis. 

In this situation, there are several cases to be distinguished: either the initial 
state of the system may be prepared according to a typical equilibrium config- 
uration, or the initial state of the system is prepared out of equilibrium, but it 
reaches some sort of equilibrium at large times. 

The first situation implies that it is known what a typical equilibrium sit- 
uation looks like. One might think that it may be approximately realized as 
follows. Choose a temperature T in the Fermi-Dirac distribution (|6.7| ) specifying 
the fermion initial conditions. Prepare the Bose fields in a classical equilibrium 
configuration at this temperature, i.e. draw them from the classical canonical 
distribution exp(— Eb/T). Correlation functions are then computed by averag- 
ing over such initial configurations. Note that in this way one has prepared two 
subsystems which are separately in equilibrium, i.e. the Bose and the fermion 
subsystem, and that these are coupled together at t > 0. If the theory is weakly 
coupled, and the back reaction of the fermions is perturbative, this can be anal- 
ysed using perturbation theory, along the lines as in 3+1 ^(J and 1 + 1 |I2"| 
dimensions for a purely classical Bose system. In the classical theory (without 
the fermions), perturbation theory leads to the conclusion [1T2J that the classi- 
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cal approximation is goodQ in 1+1 dimensions, for T = v s y/\T', with v 2 3> 1 
(cf. ( |S.1D ) and T" = 0(1). If the fermions are included, we should be able to 
reproduce numerically this way the fermion and boson propagators calculated 
analytically in the one loop approximation in the appropriate momentum and 
frequency regime. 

However, an actual numerical implementation of this idea leads us to the 
conclusion that for times beyond the validity of the perturbative calculation, the 
fermion back reaction becomes non-perturbative and the coupled system is in fact 
still far from equilibrium. 

The other possibility, the system is initialized in a nonequilibrium state, but 
reaches some sort of equilibrium at large times, relies on the assumption that the 
effective equations of motion are chaotic in the relevant coarse grained variables, 
such that at large times they lead to a distribution (by averaging over time) 
that approximates the microcanonical ensemble of the original quantum system. 
At present it is not known to us if this assumption is correct in some sense. 
We hasten to say that the objection raised sometimes to the effective equations 
that they do not incorporate scattering and therefore do not lead to the desired 
equilibrium distribution, is misleading in the inhomogeneous case considered here. 
One can easily convince oneself that there is nontrivial scattering in the effective 
equations, e.g. by taking an initial state of two fermion wave packets approaching 
each other: the localized fermions will simply experience approximately classical 
Coulomb scattering. 

Still, it is a deep fundamental question what happens at large times, in the 
original system, and in the effective field equations, given for example initially a 
bosonic configuration with a certain energy and a fermions initialized in a vacuum 
or thermal state. 

Here we only provide one numerical example to show the very long time 
behaviour. In figure [T7| we show the fermionic and bosonic energies, for initial 
Bose fields as in the previous section, i.e. ( |9.1|) , and also a similar choice of 
parameters 

N = 32, a /a = 0.1, eL = 3.2, 
v% = 8, A/e 2 = 0.25, G/e = 0.1. 

Notice, however, that the time interval < et < 45000 here is 15 times as large 
as in the previous section (0 < et < 3000). The system is evidently initially quite 
out of equilibrium, while at larger times the bosonic energy only drops very slowly. 
The bosonic energy drops by a factor of more than 15 as it gets transferred to 
the fermions. This may be related to our neglect of the bosonic fluctuations. To 
visualize the effect of decreasing bosonic energy in a typical bosonic observable, 
we show in figure [18| the charge density of the scalar field as of function of space 
x (in lattice units, x/a G {0, . . . , N — 1}) and time et. Note that initially the size 

7 In the sense of the absence of ultraviolet divergences in classical thermal loops and the 
agreement with the results from high temperature quantum field theory. 
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of the amplitudes decreases, but then remains approximately of the same size. 
Furthermore the system remains fluctuating, and inhomogeneous. More detailed 
analysis is needed to determine if the end result is approximately thermal in the 
quantum sense. 

If the system does thermalize in a reasonable way we have an approximate 
nonperturbative method for computing at finite temperature. An application 
would be e.g. the thermal sphaleron rate in the presence of CP violating fermions. 
Last, but not least, the same would apply to finite density. It is not difficult to 
prepare an initial state with a non-zero conserved quantum number (e.g. 
defined in section |3|). The corresponding finite density microcanonical ensemble 
could then emerge at large times. It would provide a new approach to the problem 
of non-perturbative finite density computations in QCD. 

11 Summary and outlook 

We considered a coupled system of classical Bose fields and a fully quantized 
fermion field, to approximate non-perturbative real-time dynamics in quantum 
field theory. The fermion back reaction was computed using a mode function ex- 
pansion, similar to previous work in the homogeneous case. However, these mode 
functions still contain the full space dependence and this makes the numerical 
problem scale like N 2d+l in d space dimensions. The immediate consequence is 
that the calculation of the back reaction is numerically much more demanding 
than in the homogeneous situation. The 1+1 dimensional model served as a 
useful laboratory for testing. 

The lattice discretization leads to fermion doubling, which is well-known in 
euclidean lattice gauge theory. We added a Wilson term in space to deal with 
the space doublers, but used 'naive fermions' in time. The time doublers were 
interpreted as a second flavour. In hindsight the treatment of fermion doubling 
could perhaps be done in a more elegant way, e.g. using staggered fermions, which 
might also lead to smaller discretization errors. 

We have studied several dynamical questions with the effective equations of 
motion. First we used external time dependent Bose fields to test the treatment 
of the fermions, using the mode functions. We found that the anomaly equation 
relating the anomalous charge and Chern-Simons number is correctly reproduced. 
However, for a finite Yukawa coupling discretization errors are larger and qual- 
itatively different in the sense that the anomalous charge looses Chern-Simons 
number after some time. This time increases with decreasing lattice distance. 

Next we studied the quantum divergence which is present in the effective 
equations of motion and we showed how this can be renormalized. We then 
demonstrated the applicability of the approximation to nonequilibrium time evo- 
lution on a few initial value problems. Finally, we have discussed the possibility 
that the effective dynamics can be used to describe a plasma at finite temperature 



30 



and density. We have shown some preliminary results at very long times which 
indicated equilibration to an ensemble of (inhomogeneous) configurations. 

Possible directions for the future are many. CP violation can be incorpo- 
rated by extending the model with more fermion fields, relevant for baryogenesis. 
Analytical calculations of various damping mechanisms of Bose fields due to a 
thermal bath of fermions can be checked numerically. In the typical examples we 
have shown, there is a large energy transfer from the Bose fields to the fermion 
fluctuations, as in inflationary scenarios. This energy transfer can presumably be 
interpreted as particle production, using the idea of adiabatic particle numbers. 
An analysis of this would require a space dependent Bogoliubov transformation. 
All of this is on reasonably firm theoretical grounds. Concerning the long time 
behaviour, and the possible approach to (some) equilibrium, we would like to 
stress that many questions are still open and require a detailed study. 

A big issue is of course how to extend this approach to 3+1 dimensions. 
Brute force using all the mode functions seems hopeless, but useful inspiration 
can perhaps be found in the dynamical fermion methods used in four dimensional 
euclidean lattice QCD. It should of course also be kept in mind that the ultraviolet 
properties of classical gauge fields in thermal equilibrium are different in 3 + 1 
dimensions than in 1 + 1 dimensions. 
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A Dirac matrices, Dirac and Majorana spinors 

Our conventions are the following. The Dirac matrices obey {7^,7"} = g^, with 
—g 00 = g 11 = 1. The antisymmetric Levi-Civita tensor e^ v is defined by eoi = 1. 
75 = — 7°7 1 . The charge conjugation matrix obeys & = C _1 ,C T = —C, and 
Cy T £t — _ an( j hence 

C 7 /C t = -75, cpI r & = P R , L , 

which is different than in 3 + 1 dimensions. We also use the hermitian matrices 
a 1 = — 7°7 1 , / 9 = 27°. The adjoint spinor is ip — ip^fi. 

As explicit representation we use the Majorana- Weyl representation 7 = 
—i&2, 7 1 = 01, 75 = 03, and C = (3. 

Working on the lattice with Wilson fermions, the following combination often 
appears 

In this paper we restrict ourselves to r = 0. 



31 



Finally, Pauli matrices acting on the 1,2 components of a Majorana spinor 
are denoted with ^1,2,3 and Pauli matrices acting in flavour space with t 12 ,3. 



B Eigenspinors in a vacuum background 

In this paper we take as initial mode functions the eigenspinors of the Dirac 
hamiltonian in a vacuum background of Bose fields, i.e. <p = v r/ V%, A\ — 0. In 
this case the eigenspinors can be found analytically, by spatial Fourier transfor- 
mation. 

The Dirac hamiltonian ( |3.17 ) in a vacuum background is given in momentum 
space by 

H D&P = a 1 s p + P(m p + p 3 T 3 m F ), (B.l) 

with s p = a -1 sinap, the spatial Wilson mass m p = a r\(l — cosap), and rriF = 
Gvr/ v2. p 3 acts on the Majorana components, and r 3 acts in flavour space. We 
use antiperiodic boundary conditions for the fermion field. This quantizes p as 

P= T {n ~\ ] ' n€{~N + l } ...^N}, (B.2) 
assuming N is even. The positive energy eigenvalues are summarized by 



Ep V = -i Is; 1 + ml , m pv = m p + rjm F , r) = ±. 



rj refers to the product of the eigenvalues of p 3 and r 3 . The mixing of the Wilson 
mass with the mass due to the Higgs vacuum expectation value gives a mass 
splitting and lifts the degeneracy. 

The eigenspinors can be written conveniently using the 2-spinor 



^2E pv (E m - Sp ) V KEm ~ s 



which is normalized to one, up v u pri = 1. 

We denote the positive energy Majorana eigenspinors, with energy E pv and 
flavours u, d, with U pr)U and U pv d- They are easily constructed, 

u P+u = ( Up + ) , u P - u =^ u °^y u P+d =(^^y u P ^ d = ^ 

The eigenspinors needed in the mode function expansion in section [3] are then 
given by 

Uauipt) 'jY^prpii U ad yX) -j=U p^. (B.3) 

V L V L 

where we used the collective label a = (p, rj). The initial conditions for the mode 
functions in section ^ are given by ( |3.21| ). 
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C Spectral flow 



The idea underlying spectral flow is that the equations of motion for the fermions 
are solved in the presence of external fields, in the special case that these exter- 
nal fields are extremely slowly varying. Then it is possible to use the adiabatic 
theorem p7| and determine the time evolution by diagonalization of the Dirac 



hamiltonian for each configuration of external fields (instead of solving the equa- 
tions explicitly). 

In this appendix we discuss the spectral flow analysis, for the special case of 
zero Yukawa coupling and a homogeneous gauge field. We use Dirac fermions, 
with the action ( |3.3| ) and G — 0. In this case the spectrum can be found analyti- 
cally for each value of the gauge field. The difference with the real-time analysis 
in the main part of this paper, is that the spectrum is calculated instantaneously 
for a given gauge field configuration. Hence the time doublers do not appear. 
The doublers in space are lifted due to the Wilson term. We use lattice units 
a = 1 in this appendix. 

Let's start with the free (A\ = 0) Dirac hamiltonian. In momentum space, it 
is given by 

H Dp = a^Sp + f3m p , 
using the notation from appendix |B[ The eigenvalues and eigenspinors are 



with 



'HDpW pe = eE p w pe , E p = yjsl + mj, e = ± 
1 



• ~ V 2E P (E P - es p ) V ^E P ~ V ' ' W 



m p 



Due to the Wilson mass term, these eigenspinors are not eigenspinors of 75 with 
eigenvalues ±1. Instead they obey 

wl e i5W P e> = e-^-5 ee >, (C.l) 

Cj p 

where s p /E p lies between -1 and +1. The sign, however, can still be used to 
indicate the would-be chirality in the continuum. 
The fermion field is expanded as0 



N 

pe 



e ipx a nf w. 



pe "J pei 



where the creation and annihilation operators obey the usual anticommutation 
relations 

,t 



{ a pe, a p ' t '} — S pp 'S e 



^Usually the w's and a's are renamed as w p + — u p , w_ p _ = v pi a p + — b p , a_ p _ = d p . 
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and zero for the other ones. The vacuum, i.e. the Dirac sea, is given by the 
non-zero expectation values 

(a p+ a P+ ) = (a p _a p _) = 1. 



In figure |19| a we show the spectrum and the filled Dirac sea for N = 8. The total 



(bare) energy and the axial charge are given by 

E = ^(^H D ^) = -Y^E P , (C.2) 



x p 

S r 



E p 

x p ^ 

The total axial charge is zero for the Dirac sea because of the precise cancella- 
tion between the + and the — branch. Because of the antiperiodic boundary 
conditions, there is no exact zero energy state. 

The case of a homogeneous gauge field is now simply found by minimal sub- 
stitution 

p= — (n — > p A] = — in 1 — C), 

y N y 2 J y 2 1 N y 2 2 h 

where we use that the Chern-Simons number C = —NAi/2tt. Recall that the 

fermions have charge q — \. In figure [19|, we show the spectrum for three values 



of the Chern-Simons number. For C — 1, the spectrum is not identical to the 
spectrum at C — 0. There are now two states with exactly zero energy (one filled 
and one empty), and there are two states exactly at the boundary of the Brillouin 
zone (p = ±7r). This is not contradictory, C + 1 is not gauge equivalent to C, 
due to the q = \ charged fermions. C + 2 is gauge equivalent to C, and indeed, 
the spectrum shown in figure p~9| c for C = 2 is identical to the original one with 
C = 0. Only the occupied and unoccupied states have been shuffled. 

The precise way in which this shuffling is done cannot be found using only 
the instantaneous spectrum. However, in the adiabatic approximation the Chern- 
Simons number has an implicit time dependence, and this gives a relation between 
the (un)occupied states for different values of C, by continuity. For instance, 
states that seem to disappear from the spectrum at (p = ±7r, E) return at (p = 
=F7r, E), because of the periodicity of the lattice. A calculation of Q§ shows that 
the contributions from the + and the — branch no longer cancel, and we find 

Q 5 (C)-Q,(0)-C-0, 

which is of course the anomaly equation. Here it is important that states that 
are close to the continuum, i.e. with \E P \ 1, do have an axial charge of 
approximately ±1. This is not the case for states deep in the Dirac sea (see Eq. 

We continue and look at larger Chern-Simons numbers. Since we only have a 
finite number of states, at some point all the occupied states have moved to the 
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+ branch. This is shown in figure pOja. In our model this happens when C = N . 
This particular distribution of occupied states has the largest axial charge that 
is possible in the lattice model, and it is given by 

Q§ma,x _ 1 \ s p\ N—*oo ^ f dp S p ^ 2 



N N ^ E p J 2vr E p vr 



where the last expression is valid in the limit N — > oo. When we continue to 
even larger Chern-Simons numbers, all the occupied states end up in the Dirac 



sky, and the Dirac sea has completely dried up, as shown in figure |20jb. The axial 
charge is zero again. 

It is clear that for these large values of the Chern-Simons number, lattice 
artefacts are substantial. 
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Figure 2: Characterization of the Bose fields configuration: the Chern-Simons 
number, the Higgs winding number and the potential energy in the Bose fields, 



normalized with the sphaleron energy, versus t', for parameters ([7.1 ). 
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Figure 3: Response of the fermions: the axial charge Q§ (solid, with rapid oscil- 
lations) and the Chern-Simons number C (dotted) versus for several sphaleron 
transitions, C = — > 4. The lines fall mostly on top of each other, in accordance 
with the anomaly equation. Parameters as in ( |7.1| ). 
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Figure 4: As in figure ||. Blow up showing C (dotted, slowly varying) and Q5 
with rapid oscillations oc Oq(— l)*/ a ° due to the time discretization, for two values 
of the lattice spacing in time, a Q /a = 0.1 (larger amplitude) and 0.05 (smaller 
amplitude). 
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Figure 5: Energy of the fermions and the sum of the energy of the fermions and 
the potential energy in the Bose fields, normalized with the sphaleron energy, 
versus t'. Parameters as in figure |3|, with a^/a = 0.005. 
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Figure 6: Lattice periodicity of Q$ for zero Yukawa coupling: Q5/N versus t'/N 
for two values of the lattice spacing, iV = 16, 32. After the rescaling with N, the 
lines converge to the same curve, as expected. Other parameters as in figure |3]. 
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Figure 7: As in figure |3], C and Q 5 for three values of the Yukawa coupling, G/e = 
0,0.1,0.5. For larger Yukawa couplings, lattice artefacts are more important. 



43 




Figure 8: As in figure [3], C and Q$ for fixed Yukawa coupling, G/e = 0.1, and 
three values of the lattice spacing, N = 16, 32, 64. Decreasing the lattice spacing 
(increasing N) shows convergence to the continuum limit. 
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Figure 9: Need for renormalization: Chern-Simons number C versus et for fixed 
bare vb, for three values of the lattice spacing, for parameters flS.4|) and initial 
conditions (|8.5[ ). The lines show dependence on the lattice spacing. 
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Figure 10: Renormalization: as in figure |9], C versus et for fixed renormalized 
vr, and three values of the lattice spacing. The bare parameter vb is determined 
from ( |8.3| ). The line with largest lattice spacing (N = 16) shows the largest 
discrepancy, and convergence towards the continuum limit is clearly visible. 
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Figure 11: Need for renormalization: as in figure the scalar field amplitude 
1 0| 2 versus et, for fixed bare vb, and three values of the lattice spacing. Again 
there is lattice spacing dependence. 



47 



7.0 



6.0 



CM 



:©: 5.0 



4.0 




Af=48,32,16 



3.0 



0.0 2.0 4.0 6.0 8.0 10.0 

et 



Figure 12: Renormalization: as in figure [Ll[ |0| versus et, for fixed renormalized 



vr, and three values of the lattice spacing. The lines show a quick convergence 
towards the continuum limit, the lines with N = 32, 48 are already indistinguish- 
able. 
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Figure 13: Non-perturbative dynamics: energy of the Bose fields Eb/e, the 
fermion field Ef/e, and the sum E tot /e versus et, for initial conditions (|9.1| ) 
and parameters (|9.2|) with t>|, = 8, G/e = 0. 
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Figure 14: As in figure [L3[ Chern-Simons number C (dotted) and axial charge Q5 
(solid) versus et. The lines cannot always be distinguished, which is in accordance 
with the anomaly equation. Small deviations, e.g., around et ~ 1300, are due to 
the finite lattice spacing. 
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Figure 15: As in figure 13, Chern-Simons number C (dotted) and Higgs winding 
number (solid) versus et. The lines are anticorrelated, which is expected when 
sphaleron transitions dominate. 
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Figure 16: As in figure [13|, Chern-Simons number C (dotted) and axial charge 
Qs (solid) versus et, for finite Yukawa coupling G/e = 0.1, and v 2 R = 4. After a 
finite time, looses C, and the anomaly equation is no longer satisfied. 
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Figure 17: Long time behaviour, energies versus et. Initial conditions and 
parameters as in figure [H| except that G/e = 0.1. 
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a) b) c) 

Figure 19: Spectrum of the Dirac hamiltonian and distribution of occupied states 
on a spatial lattice with N = 8 sites, for a) C = 0, the filled Dirac sea, b) C = 1, 
c) C — 2. The occupied states are indicated with a full dot, the empty states 
with a open dot. The response to an increasing Chern-Simons number is indicated 
with an arrow. The ±1 denotes the sign of the 75 expectation value of the states 
on the particular branch. 
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